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It was pointed out [Phys. Rev. C 91, 021302(R)] that the isotope shifts of the Pb nuclei, the 
kink at = 126 in particular, can be well described by the Hartree-Fock-Bogolyubov calculations if 
a density-dependent LS interaction derived from the 2>N interaction is incorporated. Effects of the 
density-dependence in the LS channel on the isotope shifts are extensively investigated for the Ca, Ni 
and Sn isotopes, using the semi-realistic M3Y-P6 interaction and its LS modified variant M3Y-P6a, 
as in the Pb case. It is found that almost equal charge radii between "^^Ca and '*®Ca are reproduced, 
as well as the isotope shifts in a long chain of the Sn nuclei, owing to the density-dependence in the 
LS channel. A kink is predicted at Y = 82 for the isotope shifts of the Sn nuclei, in clear contrast 
to the interactions without the density-dependence. 

PACS numbers: 21.10.Ft, 21.30.Fe, 21.60.Jz 


I. INTRODUCTION 

In the nuclear shell structure, which is formed by a 
series of the nucleonic single-particle (s.p.) orbits under 
the nuclear mean field (MF), the spin-orbit (Is) splitting 
plays an essential role. However, though known from 
the data, size of the Is splitting has been difficult to 
be accounted for only by the two-nucleon (2N) interac¬ 
tion [ij. Recent development of the chiral effective-field 
theory (yEFT) indicates that the three-nucleon (3iV) 
interaction gives significant density-dependence in an LS 
channel when it is converted to an effective 2N interac¬ 
tion ii , and that this may account for the missing part 
of the Is splitting 00. 

There have been precise measurements on isotope 
shifts 00. As well as the Is splitting, they have sup¬ 
plied problems in nuclear structure theory that have not 
been solved for decades. As an important example, the 
isotope shifts of the Pb nuclei show a conspicuous kink 
at iV = 126 when plotted as a function of the neutron 
number N. This kink was almost reproduced in a rela¬ 
tivistic MF approach 0, though not in non-relativistic 
approaches with the Skyrme interactions. This differ¬ 
ence has been recognized to originate from the isospin- 
dependence of the LS channels [9| , and has lead to exten¬ 
sion of the Skyrme energy density functional M- The 
nlggf 2 and nOiii /2 s.p. levels are relevant to the kink at 
N = 126, which are nearly degenerate in the relativistic 
model and in the non-relativistic model with modified 
LS channels of Ref. IRI. However, such degeneracy is 
not observed in ^°®Pb [^. The author, together with a 
coauthor, pointed out that the density-dependence in the 
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LS channel enables us to reproduce the isotope shifts of 
the Pb nuclei fairly well, without degeneracy between the 
711^9/2 and nOiii /2 orbitals [l^. This effect on the iso¬ 
tope shifts provides us with evidence for the 3N LS inter¬ 
action that is practically independent of the Is splitting, 
suggesting a possibility that the two problems, origin of 
the Is splitting and the isotope shifts in Pb, are simulta¬ 
neously solved by the 3iV LS interaction. In this paper 
effects of the density-dependent LS channel are exten¬ 
sively examined for the isotope shifts of other Z = magic 
nuclei. 

Whereas a density-dependent LS interaction was ar¬ 
gued in Ref. [l^, motivated to account for the model- 
dependence of the isotope shifts in the Pb nuclei, pri¬ 
mary origin of the difference between the relativistic and 
the Skyrme approaches has turned out to be the isospin- 
dependence of the LS channels. Since then the density- 
dependence in the LS channel has been considered only 
in a limited number of studies 0 , and its influence on 
physical quantities other than energies has not been ex¬ 
plored sufficiently. 


II. FRAMEWORK 

In the present work, the spherical Hartree-Fock- 
Bogolyubov (HFB) calculations are implemented for the 
Ca, Ni and Sn nuclei, as for Pb in Ref. 0. Because 
the proton number Z is magic in these nuclei, it seems 
reasonable to investigate their systematic behavior in 
the spherical HFB regime, although there could be ex¬ 
ceptions as will be mentioned below. The computa¬ 
tional method is an extensive application of the Gaus¬ 
sian expansion method [l0, and has been summarized 
in Ref. [l^. The effective Hamiltonian is comprised of 
the nuclear. Coulomb and center-of-mass (c.m.) parts: 
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H = Hn + Vc — Hc.ru.- The nuclear part is taken 
to be non-relativistic and isoscalar, comprised of the cen¬ 
tral LS tensor and the density- 

dependent central LS channels for the 

2N interaction, as well as of the kinetic energy term. In 
the former papers [iMi was denoted by 

For ii(LSp) the following form is assumed, 

^ 2i D[p(Rij)] pij X d(ry ) py ■ (s, -1- s^), (1) 

where -r^, = (r^ -hrj )/2, pij = (pi -pj)/2, 

p(r) is the isoscalar nucleon density and is the spin op¬ 
erator, with i and j being indices of constituent nucleons. 
Vc contains the exchange terms without approximation 
and Hc.m. consists both of one- and two-body terms. 

The Michigan-three-range-Yukawa (M3Y) interactions 
were obtained from the G-matrix [20|. By modify¬ 
ing the M3Y-Paris interaction [2l| phenomenologically, 
the semi-realistic M3Y-Pn interactions have been de¬ 
veloped [H, 113, whose density-independent terms 
(X = C, LS and TN) are expressed by the Yukawa func¬ 
tion. Responsible for the shell structure, the is split¬ 
ting is significant to describing nuclear structure. Ex¬ 
cept M3Y-P6a, uT®) has been enhanced from that of the 
original M3Y-Paris interaction, to reproduce s.p. level se¬ 
quence in the absence of With keeping the tensor 

channels of the M3Y-Paris interaction, Z- or N- 

dependence of the shell structure is obtained reasonably 
well, as exemplified by the level inversion of pQd ^/2 and 
plsi /2 from “^^Ca to '^®Ca [13 • Among several parameter- 
sets, it has been found that M3Y-P6 gives prediction of 
magic numbers compatible with most available experi¬ 
mental data, in wide range of the nuclear chart including 
unstable nuclei [13 • Here M3Y-P6 is taken as a yard¬ 
stick for investigating effects of the 3A^ LS interaction. 
The values of the parameters [13 in M3Y-P6 have been 
given in Ref. M- 

The yEFT derives density-dependence in the LS chan¬ 
nel as an effect of the 3N interaction Si> which may 
complement the 2N LS interaction with respect to the is 
splitting. Based on this xEFT indication, has been 

introduced in Ref. [13, instead of enhancing which 

yields a variant of the M3Y-P6 interaction called M3Y- 
P6a. However, since the currently available xEFT is not 
yet convergent at p ^ po, where po denotes the satura¬ 
tion density, the xEFT results have not completely been 
followed in quantitative respect. The functional D[p\ is 
taken to be 


The di term of the denominator, which has been fixed 
to be 1.0 fm^ in M3Y-P6a, is employed to avoid insta¬ 
bility for increasing density, whereas the results are not 
sensitive to di as will be shown in Sec. Ell Since the effec¬ 
tive interactions have been more or less adjusted to the 
observed is splitting, size of the is splitting should not 
change much even when the density-dependence is taken 


into account. Therefore the remaining parameter wi has 
been fitted to the split ting of the nOi orbits obtained 
with M3Y-P6 at ^°®Pb [ij ]. This M3Y-P6a interaction 
will be applied to other Z — magic nuclei. Note that all 
the parameters of M3Y-P6a except in a^®) and 
are identical to those of M3Y-P6. 

Since tiT®^) in Eq. CD has zero range, its contribution 
to the total energy can be represented by a functional of 
local currents. The particle-hole terms, which appear in 
the Hartree-Fock (HF) regime, are summarized as 

= i / d^^D[p{r)] {p(r) V • J(r) + ^ p,(r) V • J.(r) 

r—p,n 

-hiJ(r) •j*(r)-f-i J^(r)-j;(r) 


-iJ*(r) •j(r)-i Y JrW'jrW 

T—p^n 

-Q(r)-cr(r)- Y Qr(r) ■crr(r)j, 

T—p,n 

(3) 

where the density p(r) and the other local currents are 
defined by 


pA) = Y 

Pt ir)= Y P“/3')'pA)0a(r) , 

r—p,n 


jA) = Y 

jr(r) = -* Y Pa/3<^pA) V0a(r), 

T—p,n 


QA) = Y Q^A), 

Clr{r)=i Y Pap V0^(r) X V0a(r), 

r—p,n 

OL,(3^r 

II 

M 

(Tr{-r)=2 Y £'ap0^(r)s0a(r), 

T—p,n 


JA) = 51 J-W’ 

T—p,n 

3r{r) = 2i Y Pap 0^(r) S X V0a(r), 

Ct,0€.T 


(4) 


with the s.p. basis function 0a (r) and the one-body den¬ 
sity matrix pap = (4’|a^aa|‘&) for the HF or HFB state 
1$). It should be noticed that, because of the presence of 
D[p], integration by parts does not simplify the expres¬ 
sion of . 

Under the spherical symmetry, Qr(i’) = <7’r(r) = 
0, i[jr(r) - j*(r)] = \7pr{r) and J^(r) = j;(r) = 
Eaper PaP ^pA) (2€-s)0a(r) with £ = r X p (0a, 
0p are postulated to be spherical bases). E^^'^ yields 
the is potential as follows [^ : 


(pA)+Pr(r)) 

, , (5) 

The same Hamiltonian is applied to the pairing chan¬ 
nels in the HFB calculations, including the two-body 
term of iLc.m.- Since D[p{r)] does not contain the pair¬ 
ing tensor, contribution of to the pair energy is 

represented analogously to that given in Ref. [13 ■ 


is. 
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The semi-realistic M3Y-P6a [T^ interaction is primar¬ 
ily used in this woit and its results are compared with 
those of M3Y-P6 [l^ to exhibit effects of the density- 
dependence in the LS channel. As mentioned above, the 
strength parameter wi in M3Y-P6a has been determined 
so as to equate the nOi splitting at ^°®Pb to its counter¬ 
part in M3Y-P6. This makes influence on other is split¬ 
ting insignificant as well. For instance, the nOd splitting 
at ^°Ca decreases merely by 4 % if replacing M3Y-P6 by 
M3Y-P6a. The equal-filling approximation is applied 
to odd-A^ nuclei, whose ground states have one quasipar¬ 
ticle. It should be noted that the results of M3Y-P6 are 
qualitatively similar to those of the other interactions 
that lack density-dependence in the LS channels. To 
show it, the Gogny-DIM interaction will be taken 
as an example, which has successfully been applied to 
describing structure of many nuclei. While calculations 
beyond the spherical HFB have been implemented with 
DIM and other globally fitted interactions, qualitative 
improvement over the spherical HFB regime has not been 
reported so far on the isotope shifts of the nuclei under 
investigation, to my best knowledge. 


III. RESULTS AND DISCUSSIONS 


The frequency difference of corresponding atomic de¬ 
excitations among isotopes is converted to the difference 
in mean-square (m.s.) charge radii of nuclei. Hence the 
isotope shifts are expressed by the m.s. charge radius of 
a certain nuclide relative to that of a reference nuclide 
with equal atomic number Z. The m.s. charge radius of 
the nucleus is given by 

{r\{^Z) = {r\{^Z) + {r%{p), ( 6 ) 

where {r^)p{^Z) represents the m.s. radius for point- 
proton distribution in and {r^)c{p) the charge m.s. 
radius of a single proton. By subtracting the c.m. con¬ 
tribution, {r^)p{^Z) is obtained as 

{r^)p{^Z) = ^ ^($|(U - R)"|$) . (R = ^ E 

i^p i 

(7) 

Because {r^)c{p) is canceled out, the isotope shift is de¬ 
noted by A{r^)p{^Z) = {r^)p{^Z) — {r'^)p[^°Z) in this 
paper, with taken as a reference. 


A. Pb isotopes, revisited 

It is customary to define the isotope shift of the Pb 
nuclei by adopting ^°®Pb as a reference, Z\(r^)p(^Pb) = 
(r^)p(^Pb) — (r^)p(^°®Pb). As has been shown by com- 
paring the M3Y-P6 and M3Y-P6a results in Ref. [l^, the 
isotope shifts of the Pb nuclei can be improved with the 
density-dependent LS interaction. 
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FIG. 1. (Color online) Isotope shifts of the Pb nuclei 
A(r^)p('^Pb) obtained from the HFB calculations, (a) Com¬ 
parison among the results with M3Y-P6 (green dashed line), 
M3Y-P6a (red solid line) and its variant of di = 2.0 fm® 
(black dot-dashed line), (b) Comparison between the resnlts 
with DIM (brown dot-dashed line) and DlMa (orange short- 
dashed line). Experimental data qnoted from Ref. [sj] are 
shown by crosses for reference. 


It is exhibited in Fig. [TJa) that the di parameter in 
Eq. ([2]), which has been presumed to be 1.0fm^ in M3Y- 
P6a, hardly influences the results. If di is assumed to 
be 2.0 fm^, wi = 842 MeV • fm® follows by equating the 
nOi splitting to that with M3Y-P6 at ^°®Pb, instead of 
742MeV • fm® in M3Y-P6a. Notice that D[po\ values are 
close to each other. The results with di = 2.0 fm® are 
shown by the black dot-dashed line. Comparison to the 
results with M3Y-P6a confirms insensitivity to di. 

The density-dependent LS term in the effective inter¬ 
action may produce different behavior in the isotope shift 
unless it is too weak, as typically represented by the kink 
at V = 126 in A{r^)p{^Vh). This is not constrained to 
the M3Y-P6 case, although it depends on the interaction 
to certain degree how strong its effects are. It is here illus¬ 
trated by DIM. The LS channel of the DIM interaction 
is obtained if D[p\ in is replaced by a constant wq. 

The strength of (z.e. wq) was determined in a fully 
phenomenological manner. When is introduced, 

it is reasonable to reduce wq so as to keep size of the 
is splitting, which leaves the ratio wipo/wo almost arbi- 
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FIG. 2. (Color online) Isotope shifts of the Ca nuclei 
Zl(r^)p(^Ca), obtained from the HFB calculations with M3Y- 
P6a (red solid line), in comparison to those with M3Y-P6 
(green dashed line) and DIM (thin brown dot-dashed line). 
Experimental data are quoted from Ref. [3i[ (crosses). 


trary. To see qualitative effects of wg is reduced 

to 30 % of its original value and wi is adjusted by equat¬ 
ing the nOi splitting to that of the original DIM, with 
taking di = 1.0fm^. This interaction is called DlMa in 
the present article. The results of DlMa on the isotope 
shifts of the Pb nuclei are compared to the DIM results 
in Fig. [ijb). The experimental data are shown for refer¬ 
ence, although comparison to them is not the point here. 
As has been pointed out [13,[SI, the slope of the isotope 
shifts becomes steeper in > 126 than in < 126 be¬ 
cause of the occupation on nOiii/ 2 - With DIM (DlMa), 
nlgg /2 lies lower than nOzii /2 by 1.4 MeV (1.2 MeV) in 
208pb energy difference suppresses occupation on 

nOjii /2 near ^°®Pb, and therefore the effect of the density- 
dependent LS interaction is hindered in N > 126, though 
present. However, because the nOiis/g function shrinks, 
the isotope shift varies more slowly with DlMa than with 
DIM in A^ < 126, producing a visible kink at A^ = 126. 


A(r^}p(^^Ca). As illustrated for nOi orbits in Fig. 1 
of Ref. [13 , the density-dependence in the LS channel 
tends to shrink the j = £ -I- 1/2 orbits while extends the 
j = £ — 1/2 orbits. Therefore the m.s. radius of neu¬ 
trons occupying the O/ 7/2 orbit comes smaller, decreas¬ 
ing by 0.23 fm^ at ^°Ca when M3Y-P6 is switched to 
M3Y-P6a. Moreover, the m.s. radius of nO/7/2 increases 
from 40Ca to ‘‘SCa only by 0.035 W with M3Y-P6a, in 
contrast to 0.365 fm^ with M3Y-P6. This is interpreted 
as the neutrons on O/ 7/2 feel mutual attraction and the 
density-dependent LS channel blocks them to distribute 
more broadly, which further suppresses These ef¬ 

fects lead to the result Z\(r^)p(‘^®Ca) ss 0. Thus it may 
be possible to ascribe the very small difference in the 
charge radii between "^^Ca and "'^®Ca to an effect of the 
3A^ interaction. 

Though not shown to keep the figure visible, 
the density-dependence in the LS channel makes 
Z\(r^)p(‘^®Ca) smaller also when switching DIM to 
DlMa, confirming the qualitative effect of the 
density-dependence. In quantitative respect, however, 
Z\(r^)p(‘^®Ca) with DlMa is close to that with M3Y- 
P 6 , mainly because DIM gives almost twice larger 
Z\(r^)p("^®Ca) than M3Y-P6. 

The isotope shifts of show sizable deviation 

from those of Since (r^)p("^Ca) varies almost 

linearly from “^^Ca to "^^Ca in the spherical HFB calcu¬ 
lations, none of the interactions including M3Y-P6a re¬ 
produce the observed A^-dependence of Z\(r^)p(^Ca) in 
N = 22 — 26. This discrepancy may be ascribed to ef¬ 
fects beyond MF, including influence of the a-clustering. 
On the contrary, the rising of A{r^)p{^Ga) from ^®Ca to 
®°Ca in the M3Y-P6a result is in good agreement with the 
data, while not so good in the M3Y-P6 and DIM results. 
This also supports presence of the density-dependent LS 
interaction. 


B. Ca isotopes 

The isotope shifts of the Ca {i.e. Z = 20) nu¬ 
clei are defined by taking ^°Ca as a reference nuclide, 
Z\(r2)p(^Ca) = (r2)p(^Ca) - (r2)p(40Ca). The HFB re¬ 
sults are depicted and compared with the data in Fig. [2] 
It is commented that the weak instabilities against the 
octupole [Hill and the pairing [ 2 I correlations in “^^Ca 
are lifted with M3Y-P6a. The inversion of the pOd ^/2 
and plsi /2 levels from "‘“Ca to "‘^Ca is reproduced with 
M3Y-P6a as well as with the previous interactions [^ . 
in which plays a significant role. 

As doubly magic nuclei, both "‘“Ca and ^®Ca are ex¬ 
pected to be well described within the HF regime. Al¬ 
though it has been known that their charge radii are close 
to each other, this property has been difficult to be re¬ 
produced by self-consistent nuclear structure calculations 
so far. It is interesting to see A(r^)p("‘®Ca) ~ 0, i.e. 
(r^)p("‘“Ca) ss (r^)p(^®Ca), in the M3Y-P6a result. At¬ 
traction from the neutrons occupying O/ 7/2 determines 


C. Ni isotopes 

For the Ni (i.e. Z = 28) isotopes ““Ni is used as 
a reference, defining isotope shifts as zl{r^)p("^Ni) = 
(r^)p(^Ni) — (r^)p(““Ni). The results are depicted in 
Fig. [31 

Steeper slope is obtained for A(r^)p(""‘Ni) with M3Y- 
P 6 a than with M3Y-P6 in 28 < iV < 40. This 
is attributed to contribution of neutrons occupying 
nO/5/2, which distribute more broadly when the density- 
dependence is incorporated in the LS channel. Compared 
to M3Y-P6, M3Y-P6a gives slightly better agreement 
with the measured value for “®Ni while slightly worse for 
“^Ni. Since both of the M3Y-P6 and M3Y-P6a results are 
in reasonable agreement with the existing data, it is not 
obvious to tell whether the density-dependent LS channel 
improves the theoretical results on Z\(r^)p("^Ni) or not. 
It is desired to measure neutron-rich or neutron-deficient 
isotopes, e.g. ““Ni. 









FIG. 3. (Color online) Isotope shifts of the Ni nuclei 
Z\(r^)p(^Ni). See Fig. [2]for conventions. 


TABLE I. Root-m.s. charge radii of the reference nuclei. 
Spherical HFB results with M3Y-P6 and M3Y-P6a are com¬ 
pared with experimental data [^. (/rp) and (ep) correspond^ 
to the data on (r^)c(p). 


nuclide 


M3Y-P6 

M3Y-P6a 

Exp. 

160 

(mp) 

2.732 

2.725 

2.6991T0.0052 


(ep) 

2.743 

2.737 


«Ga 

(mp) 

3.491 

3.487 

3.4776±0.0019 


(ep) 

3.500 

3.496 


®°Ni 

(PP) 

3.825 

3.779 

3.8118±0.0016 


(ep) 

3.834 

3.787 



(PP) 

4.653 

4.627 

4.6519±0.0021 


(ep) 

4.659 

4.634 


208 Pb 

(PP) 

5.493 

5.464 

5.5012±0.0013 


(ep) 

5.499 

5.470 




N 


N < 82, an effect of the density-dependent LS channel on 
n0hii^2i makes the kink more conspicuous. This kink in 
Z\(r^)p(^Sn) at = 82 is a clear contrast to the results 
of M3Y-P6 and DIM which have no density-dependent 
LS channel. Thus the density-dependence in the LS chan¬ 
nel is essential to the kink. Future measurements with 
respect to this kink {i.e. Zi(r^)p(^Sn) beyond N = 82) 
will be intriguing, which could be a touchstone of the 3N 
LS effects. 


E. Absolute values of charge radii and 
neutron-skin thickness 


FIG. 4. (Color online) Isotope shifts of the Sn nuclei 
Z\(r^)p(^Sn). See Fig. [2]for conventions. 


D. Sn isotopes 

For the Sn (i.e. Z = 50) isotopes ^^°Sn is adopted as 
a reference nuclide, by which isotope shifts are defined 
as Z\(r^)p(^Sn) = (r^)p(^Sn) — (r^)p(^^°Sn). The results 
are presented in Fig. U 

It is found that Z\(r^)p("^Sn) is well described with 
M3Y-P6a in a long chain of the Sn isotopes, owing 
to the density-dependence in the LS channel. The 
density-dependence gives narrower distribution of neu¬ 
trons occupying the Ohii /2 orbit, whose attraction re¬ 
duces Zi(r^)p('^Sn) in TV > 70. The same trend is ob¬ 
tained with DlMa. 

It is remarked that a kink is predicted at TV = 82 
with the M3Y-P6a interaction. Origin of this kink is ac¬ 
counted for analogously to the Pb case. Attraction from 
neutrons occupying 0/ig/g contributes to the broader pro¬ 
ton distribution. While occupation probability on nOhg /g 
is negligibly small in < 82, nO/ig /2 is partially oc¬ 

cupied in A > 82 owing to the pairing correlation, al¬ 
though the lowest s.p. level above A^ = 82 is nl/ 7 / 2 . 
The relatively large m.s. radius of nOhg/ 2 , for which the 
density-dependence in the LS channel is responsible, pro¬ 
duces the kink at N = 82. The reduction of the radii in 


Although the isotope shifts are the main subject of this 
paper, it deserves discussing absolute values of the charge 
radii of the reference nuclei as well. While the finite- 
size effects of protons are canceled out in the isotope 
shifts, absolute values of the charge radii are affected by 
the ambiguity in the c harg e radius of an isolated proton. 
A recent compilation [33| gives two values, {r‘^)c{p) = 


0.84087 ± 0.00039 by pp Lamb shift and 0.8775 ± 0.0051 
by ep CODATA. The spherical HFB results on the root- 
m.s. charge radii \/ (r^)ci^Z) with M3Y-P6 and M3Y- 
P 6 a are tabulated in Table HI obtained from Eq. (© by 
taking both of the pp and the ep values for {r'^)c{p), in 
comparison to the experimental data M- 

It is found via comparison of the M3Y-P6a results to 
the M3Y-P6 ones that the density-dependence in the LS 
channel tends to reduce the radii. Whereas this effect 
is weak in and ^°Ca which are so-called LS-closed 
nuclei, it is stronger in ®°Ni, ^^°Sn and ^°®Pb, with 
0.02 — 0.05 fm difference between the two interactions. It 
is reasonable that the effect is minimal in the LS-closed 
nuclei in which all the is partners (the j = i zL 1/2 or¬ 
bitals) are filled, in contrast to the jj-closed nuclei in 
which there is a pair of an occupied j = I 1/2 and 
unoccupied j = i — 1/2 orbitals. Since the parameters 
of the central channels have been more or less fitted to 
the measured radius of ^°®Pb, the charge radius is well re¬ 
produced with M3Y-P6 for this nucleus, and so for ^^°Sn. 
This indicates that M3Y-P6a slightly underestimates the 
charge radii of ^^°Sn and ^°®Pb, although it gives values 
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closer to the measured ones than M3Y-P6 for and 
■^^Ca. As the MF approaches are expected to be the more 
appropriate for the heavier nuclei, there may be a room 
to readjust the parameters in the central channels. Such 
readjustment could influence the saturation density by 
1 %. 

The neutron-skin thickness, which is usually repre¬ 
sented by the difference of the neutron and proton root- 
m.s. radii ^— \J{r‘^)p in a single nucleus, attracts 
interests because it correlates to variation of the sym¬ 
metry energy for increasing density [s^. Because of 
their neutron excess and of experimental accessibility, 
the neutron-skin thickness is frequently argued for '*®Ca, 
®®Ni, ^^^Sn and ^°®Pb. The density-dependence in the 
LS channel hardly influences this quantity at ^°®Pb and 
at ^^^Sn; difference between M3Y-P6 and M3Y-P6a is 
0.001 fm or less. On the other hand, \J{r^)n — 
is larger by 0.011 fm with M3Y-P6a than with M3Y-P6 
at ®®Ni and smaller by 0.014 fm at "‘^Ca. This is be¬ 
cause the effects of the density-dependence tend to be 
canceled if both of the £s partners are occupied. Though 
not very strong, influence of the density-dependence in 
the LS channel cannot be discarded for certain nuclides 
in order to argue their neutron-skin thickness to 0.01 fm 
accuracy. 

IV. SUMMARY 

Effects of the density-dependent LS interaction, which 
was suggested from the chiral 37V interaction, have been 
extensively studied on the isotope shifts of the Ca, Ni 
and Sn nuclei. Since the density-dependence makes the 
LS interaction stronger in the nuclear interior relative to 
its exterior, the s.p. wave functions of the j = £ -|- 1/2 
orbits shrink while those of the j = ^ —1/2 orbits become 
broader, if the is splitting is maintained. It was pointed 
out in the previous paper that, owing to this mechanism, 
the kink in the isotope shifts of the Pb nuclei a.t N = 
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